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We consider the two-dimensional A/" = (2, 2) Wess-Zumino model with a cubic superpotential at weak and 
intermediate couplings. Refined algorithms allow for the extraction of reliable masses in a region where per- 
turbation theory no longer applies. We scrutinize the Nicolai improvement program which is supposed to guar- 
antee lattice supersymmetry and compare the results for ordinary and non-standard Wilson fermions with those 
for SLAC derivatives. It turns out that this improvement completely fails to enhance simulations for Wilson 
fermions and only leads to better results for SLAC fermions. Furthermore, even without improvement terms 
the models with all three fermion species reproduce the correct values for the fermion masses in the continuum 
limit. 
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I. INTRODUCTION 

Supersymmetric models have drawn much attention over the 
past decades. In particular, supersymmetric extensions of the 
standard model have become a primary research topic for 
model building. The additional symmetry of these models 
proves to be a very useful tool for the study of their pertur- 
bative and non-perturbative aspects. It is notoriously compli- 
cated to check and extend the predictions made by supersym- 
metry in a strong coupling regime where standard perturbation 
no longer applies. 

At the same time, lattice simulations of quantum field the- 
ories have been very successful in an increasing number of 
applications. In some theories, it is possible to match numeric 
results at weak coupling to perturbative continuum results; at 
stronger coupling, lattice simulations are often the only vi- 
able way to investigate non-perturbative properties of the the- 
ories. As nonperturbative effects are automatically taken into 
account, it is desirable to apply the lattice approach also to 
supersymmetric theories. This has been the subject of a num- 
ber of publications, see, e.g., j^ |2l and references therein. 
There are a number of challenges with respect to this goal 
since it is well known that full supersymmetry can not be re- 
alized in a generic lattice model. The reason for this can be 
traced back to the failure of the Leibniz rule on the lattice yl. 
The full supersymmetry can only be recovered in the limit of 
vanishing lattice spacing (continuum limit); but only in some 
cases, the conditions for such a restoration are under control. 
E.g., it has been shown that even in supersymmetric quantum 
mechanics the naive discretization does not lead to a super- 
symmetric continuum limit 101; generically, such a limit can 
at best be achieved by finetuning the bare coefficients of all 
supersymmetry-breakingcounterterms [50. This, however, re- 
quires much knowledge of the theory in advance. In some 
cases the relevant operators can be determined perturbatively, 
cf. B. A possible way beyond perturbation theory is the 
application of a blocking transfromation as in ^ for a free 
theory. This may lead to a solution similar to the Ginsparg- 
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Wilson relation for the chiral symmetry ||8[] . 

It is possible to reduce the number of relevant operators in 
the continuum limit if some symmetries of the continuum the- 
ory are already realized in the lattice action. The standard lore 
is that it is sufficient to realize just a part of the supersym- 
metry on the lattice in order to ensure the correct continuum 
limit. There have been many suggestions and numerical in- 
vestigations with respect to such a partial realization of the 
supersymmetry algebra on the lattice, e.g. 121 and lllOll . An 
elegant suggestion uses a Nicolai map lUlll to create lattice 
improvement terms t hat g uarantee a partial realization of su- 
persymmetry, cf. e.g. [12(1. 

Numerical simulations of supersymmetric theories face the 
further difficulty that bosons and fermions on the lattice 
should be treated on equal footing. This demands for dy- 
namical fermions; however, such simulations are notoriously 
numerically involved. Therefore, it is advisable to start with 
low-dimensional theories in order to gain information about 
the performance of the different supersymmetric lattice for- 
mulations. On the other hand, such dynamical fermion sim- 
ulations in low dimensions are interesting in their own right 
because they allow for an explicit investigation and improve- 
ment of the corresponding known algorithms. 

We have started the a naly sis of such low-dimensional mod- 
els in a previous paper lll3ll with investigations of various lat- 
tice formulations of supersymmetric quantum mechanics and 
first tests of the two-dimensional Wess-Zumino model at weak 
coupling. Here we will extend the analysis of the latter theory 
using far more elaborate numerical techniques to reach inter- 
mediate to strong values of the coupling. We are able to sim- 
ulate the Wess-Zumino model for a much larger parameter re- 
gion as in related previous works II12I1 and 1114)]. Starting from 
the standard hybrid Monte Carlo algorithm II15I1 we employ 
a novel combination of algorithms involving both a higher- 
order lllql integration scheme and Fourier acceleration lUTIl . 
This entails much better statistics in combination with larger 
lattice sizes. These improvements lead to reliable new results 
even at stronger coupling where considerable deviations from 
perturbative predictions, e.g., for the masses of the supersym- 
metric partners can be observed. 

A further goal was a systematic study of the effects of the 
above-mentioned improvement terms introduced by the Nico- 
lai map 11211 . In this paper, we present the first explicit com- 



parison of the models with and without such terms. It may 
come as a surprise that for Wilson fermions the "improvement 
term" even fails to improve the properties of the lattice model. 
Moreover, such terms introduce new complications and can 
lead to unreliable numerical results. 

In previous works lll3[ [isll it has been demonstrated that 
lattice models based on the SLAC derivative lll9|] and on the 
twisted Wilson formulation (as introduced in II13I1 ) are par- 
ticularly well-behaved as far as the continuum limit is con- 
cerned. Even at large lattice spacing the continuum result is 
approximated very well. In the current simulation the SLAC 
derivative again proves to be the best choice because it allows 
for much larger values of the coupling constant, and only a 
comparably coarse lattice is needed to extract the correct con- 
tinuum results. It is interesting to note that contrary to a real- 
ization with Wilson fermions the improvement terms for the 
SLAC derivative in fact lead to better numerical results. 

The paper is organized as follows: We start out with a 
short introduction of the different lattice realizations of the 
two-dimensional J\f = 2 Wess-Zumino model and the corre- 
sponding improvement terms with their respective lattice and 
continuum symmetries. Then, we present the numerical re- 
sults of our simulations; in particular, we compare the masses 
of the supersymmetric partners as a measure for how well su- 
persymmetry is realized on the lattice. A comparison of the 
various models with the perturbative continuum prediction at 
smaller values of the dimensionless coupling is the subject of 
Section Unci At last, we turn special attention to the regime 
of intermediate couplings where the measured masses differ 
considerably from the one-loop results. 



II. LATTICE MODELS 

A. Supersymmetrically improved lattice actions 

The lattice models under consideration have been discussed at 
length in llOl . Therefore, we shall only briefly recall the def- 
initions of the corresponding lattice actions. In terms of com- 
plex coordinates z and z for the two-dimensional Euclidean 
spacetime together with the corresponding holomorphic and 
anti-holomorphic differentials d and d the continuum action 
of the Af = 2 Wess-Zumino model reads 



'-'cnnt 



d^x{2d'fdip+^\W'{ip)\^+^ljMi(j). (1) 



The bosonic potential is given by the absolute square of the 
derivative of the holomorphic superpotential W{ip) w.r.t. its 
argument ip = ipi+iLp2- Apart from the standard kinetic term 
for the (two-component) Dirac spinors, the Dirac operator M 
contains a Yukawa coupling. 



Taken together they satisfy the A^ = (2, 2) superalgebra, and 
it has been argued that at most one supersymmetry can be 
preserved on the lattice lll2ll . With the help of the explic- 
itly known form of the Nicolai map it is possible to construct 
such a lattice model straightforwardly. In terms of the Nico- 
lai variable (^^ = 2{d(p)x + Wx on the lattice, the discretized 
Wess-Zumino action reads 



•5=5^ ixS.X + ^ ll^xMxyll^y. 



(3) 



x,V 



Here, Wx is taken to be the lattice counterpart of the contin- 
uum operator W'{ip), i.e. Wx — W'{lPx)- The matrix A/ is 
given by 



Mxy 



VV xy ^C/xy 



9Vy 



dipy 



2d W I \ -^^ -^^ 



(4) 



We require all lattice difference operators to be antisymmetric, 
dxy = —dyx- From the second equality in (HI we can read off 



that W, 



xy ■' 



dWx/dipy 



One easily checks that (O is invariant under the following 
(supersymmetry) variation. 



Sfx 

5(px 






5i!2,x = 



' 2>iX^-) 



Hl.x = 0, 
H2.X = 0. 



(5a) 
(5b) 



In terms of the original fields, (O takes the form 



S- 



=E(2(^^ 



{d'p)x + \\Wx\ 



Wx{d>f)x + Wx{d^y. 



+ 



2_^{'^l,x,ll}2,x) 



*^ xy ^C/xy 
2dxy Wxy 



i'i,y 

^2..y 



(6) 



This supersymmetrically improved lattice action differs from 
a straightforward discretization of ^ by 



A^ = ^ (Wx{dp)x + Wx{d^)x 



(7) 



a discretization of a surface term in the continuum theory 
(which is therefore expected to vanish in the continuum limit 
for suitably chosen boundary conditions). For the free theory 
{Wx = rnifx) AS* — readily follows from the antisymmetry 
of the difference operator dxy while for interacting theories 
(|2ll guarantees the invariance of the action under (|5]l without 
the need of the Leibniz rule. To study the impact of SUSY 
improvement we will compare also the improved action with 
the unimproved straightforward discretization of ([1]) (without 
AS*). 



M = Yd + Yd + W"P+ + W P. 



(2) 



In (in we have introduced chiral projectors P± = ^(1 ± 73) 
whichintheWeyl basis with 7^ = ai,Y = ^o'2, 73 = 17^7^ 
project onto the upper and lower components of ip. In the 
form ([T]) the action is invariant under four real supercharges. 



B. Lattice fermions 



For the symmetric difference operator 



(8) 



doublers are inevitably introduced into both the bosonic and 
fermionic sector In order to get rid of them a Wilson term 
may be added to the superpotential so as to maintain the in- 
variance of the action under (|5]). Within this context two dif- 
ferent choices have been discussed previously 11 1311 . 



and 



W^ = W'{^^) - -(Ay.), 



W^ = W'{ip.,) + -(A^), 



(9) 



(10) 



We stress that for Wilson fermions, the derivative of the super- 
potential is now shifted as compared to the situation after (O. 
From the first expression we recover the standard Wilson term 
for the fermions, i.e. Wxy ~ W"{ipx)5xy — ^Axy The opera- 
tor Axy is the usual two-dimensional (lattice) Laplacian 2dd. 
The second possibility ( fTOl i leads to Wxy = W" {(px)Sxy + 
j3^Axy Here, the appearance of 73 motivates the name 
twisted Wilson fermions (not be confused with the recently 
introduced twisted mass formulation of lattice QCD). It was 
already shown for the free theory foy that twisted Wilson 
fermions suffer far less from lattice artifacts than their stan- 
dard Wilson cousins. Here we will show that they remain su- 
perior even for (strongly) interacting theories. 

Besides these two (ultra-)local difference operators we have 
previously suggested to reconsider the non-local SLAC lattice 
derivative in the context of lattice Wess-Zumino models. The 
matrix elements of the SLAC derivative are most conveniently 
given for a one-dimensional lattice with an odd number of 
lattice points L, 



dx 



ity 



{-if-y- 



■k/L 



sin(7r(x — y)/L) ' 



dx 



0. 



(11) 



The generalization to higher dimensions is straightforward 
and amounts to forming suitable tensor products of (fTTl i.' For 
SLAC fermions no further modifications to the superpoten- 
tial are necessary. It is due to this fact that they constitute an 
interesting alternative to Wilson fermions. 



C. Discrete symmetries 

For the numerical analysis of Sec.|lII]we have chosen the su- 
perpotential 



W{ip) 



\mip'^ 



hv^ 



(12) 




no. 1: Classical bosonic potential V{ip) = ^\W'{ip)\'^ from OU 
shown for vanishing imaginary part (<^2 =0). In the free theory 
limit {g -^ 0) the left minimum is pushed towards minus infinity. 



allows for discrete symmetries Zf x Zj which act as reflec- 
tions interchanging the two vacua and as complex conjuga- 
tions on the complex scalar field: 



Z? 



2- 'P' 



m 
9 



^ 



and 



ip^ip. 



(13) 



so that also the potential ^|Vl^'(<y5)p is invariant under both 
transformations, cf. Fig.[T] 

From the explicit form of the fermion matrix M and its 
adjoint M^ 



M = 


-i^'^^, + m + 25(<pi + 273(^2), 


(14a) 


A/t = 


-7^5^ + m + 2.g((y9i - 173932) 


(14b) 


one finds that 






Z^ il/ H^ 


~ 73 Af 73 , Z^ : M ^ 73^/^73 , 


(15) 



which shows the invariance of the determinant.^ 

Apart from Lorentz transformation, the continuum model is 
(irrespectively of the concrete form of the superpotential) also 
invariant under time reversal and parity transformations 



Zg: {z,z) 1-^ (-Z, -z), 



Zl: 



{z,z) ^ {z,z). 



(16) 



Barring possible Wilson terms, the unimproved lattice models 
obviously inherit all discrete symmetries from the continuum. 
By contrast, the supersymmetrically improved lattice models 
are invariant only under a combination of all symmetries. We 
find 



which coincides with that in earlier simulations of the Wess- 
Zumino model III2I UM . We will assume the coupling con- 
stants 771 and g to be real and positive. The superpotential (fT2] | 






A 



-"^'xiPv).. 
W'x{d^)x. 



(17a) 
(17b) 



' The reason for an odd number of lattice points originates from a reality 
condition on the matrix elements lilt . As such it is a mere technicality in 
order to ease numerical simulations. 



■ This is true at least up to an irrelevant sign. On the lattice the fermion 
matrix M always has an even number of rows and columns, hence this 
phase does not appear. 



TABLE I: Comparison of various lattice models w.r.t. their symmetries. All statements refer to to the interacting theory, i.e. g 7^ 0. The notion 
iF^ denotes the combined action of a field and parity transformation as discussed in the text. 





(1) 


(2) 


(3) 


(4) 




(5) 




Wilson impr. 


Wilson unimpr. 


twisted Wilson" 


SLAC impr. 




SLAC unimpr. 


lattice derivative 


local 


local 


local 


non-local 




non-local 


lattice artifacts 


0{a) 


0{a) 


O^af 


'perfect' 




'perfect"^ 


modifications to superpot. 


yes 


yes 


yes 


no 




no 


discrete symmetries 


Zf 


Z^ X Z^ X 1% 


zf 


ZF^ X TF^ 


TSi 


X Z^ X Z§ X Zf 


supersymmetries 


one 


none 


one 


one 




none 



"Only improved considered. 

''In the interacting case the good scaling properties are lost. However the 
overall size of lattice artifacts is still much smaller when compared to Wilson 
fermions. 

'The dispersion relation is up to the cut-off the same as in the continuum. 



Thus, for the improved models (with SLAC fermions) the con- 
tinuum symmetry is reduced. 



Z2 X Z2 X Z2 X Z2 



sr 



xZl 



PC 



diag(Z| X Z^ X Z^) x diag(Z^ x Z^). (18) 



PC 



Here, the diagonal subgroup diag(Z2 x Z2) is a group Z2 
generated by the product of the generators of Zj and Zj (anal- 
ogous notations are used for the other groups). It readily fol- 
lows that the improvement term must have a vanishing expec- 
tation value in the original ensemble without improvement. 
We have checked this with a large numerical precision. For 
Wilson and twisted Wilson fermions with improvement the 
r. h. s. of ( fTSl l is even further broken down due to the presence 
of the (twisted) Wilson term in the superpotential. For Wilson 
fermions, the bosonic action can be read off from (|6]l and Q, 



^B = ^l]|(5^)x + M^^-i(A^). 



(19) 



Since A^^^ is invariant under both time reversal and parity, 
([TTT i cannot be preserved; the Wilson term inevitably changes 
sign. Conversely, from the bosonic action with twisted Wilson 
fermions 



^B=iEk^^)-+^^+f(^^)- 



(20) 



only (lp —^ —m/g — (p, d ^ —d) can be shown to yield a 
symmetry. In either case the breaking of the other symme- 
tries is induced by a higher-dimensional operator and may be 
expected to be at most 0{a) III2II20I1 . Nevertheless, at finite 
lattice spacing, the physics might be affected since the overall 
size of the breaking terms is a dynamical question. By con- 
trast, SLAC fermions with the larger symmetry ( [TtI i are again 
favored. 

In Tab. U we summarize all lattice models to be dealt with 
in the next section. 



III. NUMERICAL RESULTS 

As outlined in the introduction we have employed the stan- 
dard hybrid Monte Carlo algorithm for our numerical sim- 
ulations. The fermion determinant was estimated stochasti- 
cally utilizing real pseudo-fermion fields. The reason for real 
pseudo-fermions derives from the presence of only a single 
flavor such that the square root of the pseudofermionic kernel 
Q^^ = {MM'^)^^ is actually needed. We note in passing 
that the pseudo fermion action remains real with this choice 
since also the fermion matrix is real for Majorana basis is cho- 
sen. Hence the latter was adopted for all our simulations. A 
significant gain was achieved by combining higher order in- 
tegrators with Fourier acceleration techniques. With the help 
of the former one can avoid the requirement for ever smaller 
time-step sizes during the MD step of the HMC while a care- 
ful tuning of the latter allows for autocorrelation times r < 5 
over the whole range of parameters analysed. In particular for 
small lattice spacings, i.e. at large lattice sizes this was seen 
to reduce significantly critical slowing down as also reported 
in 1I21I1 . A detailed account of the algorithm employed here 
will be published separately at a later time. 



A. Dynamical properties of improved lattice actions 

Before discussing measurements of physical observables in 
the next section we will first focus on the improvement 
term (|7]|. The aim is to understand the difference between 
improved and unimproved lattice models w.rt. predictions of 
supersymmetry. One possible test is a measurement of the 
bosonic action itself. With the help of the Nicolai map ap- 
pearing in (O one can show that 



(^b) = N. 



(21) 



Here, N = N^ x N^ denotes the total number of lattice points, 
and ( I2TI 1 is only expected to hold when fermions are included 
dynamically. Then, however, this prediction holds irrespec- 
tively of the concrete value of the coupling constants. With a 
slightly different argument the same was also found in lUT 
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FIG. 2: Normalized bosonic action as a function of the bare 
mass lattice parameters using Wilson fermions with the improved 
(filled squares) and unimproved (empty squares) actions from ei- 
ther quenched (red) or dynamical fermion (blue) simulations (A'^ — 
16 X 16). 



1.08 


- 














1.06 


- 












- 


1.04 


- 












- 


1.02 


- 












- 






















A 


ft 


A 








0.98 


- 




ft 


A 




A 


A- 


0.96 


- 






A 


A 
A 


A 


A" 


0.94 


- 










A 


A~ 


0.92 
n on 


- 


1 


1 


1 


1 


1 


1 



0.0 



0.1 



0.2 



0.3 



0.4 



0.5 



0.6 



FIG. 3: Same as in Fig.|2]but for SLAC fermions. 

Equation (l2Tl i provides a test observable distinguishing im- 
proved from unimproved lattice models as well as quenched 
from dynamical fermion simulations. To accomplish this, we 
have run simulations with both (standard) Wilson and SLAC 
fermions. The results are shown as a function of the bare lat- 
tice mass parameter mun = rn/N^. Since the continuum limit 
for this theory is obtained from m\^a -^ 0, smaller values 
of miatt Ukewise mean a finer lattice spacing (and for fixed 
N a smaller spacetime volume). The dimensionless coupling 
strength A — g/m was set to A — 1. The lattice sizes 
we used for our numerical simulations were A^ = 16 x 16 
for Wilson and A^ = 15 x 15 for SLAC fermions. For the 
quenched simulations 500,000 (independent) configurations 
were evaluated, and 30,000 configurations with dynamical 
fermions were analysed. The results are shown in Fig. |2]for 
Wilson and in Fig. |3]for SLAC fermions. One clearly ob- 
serves that the quenched data significantly deviate from the 
predicted value which illustrates the necessity of dynamical 
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FIG. 4: MC history of the lattice mean $i = A''"^ J^^ ipi^x and size 
of the improvement term for Wilson fermions (A'^ = 16 x 16, A = 
0.6,miait = 0.3). 



fermion contributions in order to retain supersymmetry. Us- 
ing an unimproved action with dynamical fermions we find 
much smaller deviations which in case of the Wilson fermions 
are already hard to distinguish from the improved results. For 
SLAC fermions the deviations are somewhat more systematic 
and remain also clearly distinguishable from other dynami- 
cal fermion simulations. A second difference between Wilson 
and SLAC fermions may be infered from Fig. HI Namely, 
there is a distinct correlation between the ground state around 
which the field (pi fluctuates on the one hand and size and 
variance of the improvement term on the other hand. This 
may be taken as direct manifestation of the additionally bro- 
ken Z|^'^-symmetry due to the Wilson term which will also 
play a role when discussing the failure of improvement in the 
next paragraphs. 



Limitations of improved lattice actions 

Studying the improvement term AS for models with either 
Wilson or SLAC fermions we have found that the system is 
ultimately pushed into an unphysical region of configuration 
space, at least for strong coupling. Our simulations have re- 
vealed that this instability is controlled by the actual size of 
the bare mass parameter and the coupling strength A. Simu- 
lations tend to fail more often as either of them grows. The 
study of this phenomenon with Wilson fermions turns out to 
be clumsy since there is no clear correlation between the value 
of the coupling and the number of configurations where the 
instability occurs. Hence we prefer to present our analysis 
from the simulations with SLAC fermions. However, it should 
be emphasized again that for either Wilson or twisted Wilson 
fermions the qualitative picture is the same as described be- 
low. 

It is to be expected that the improvement term grows with 
the coupling strength A and vanishes continuously in the con- 
tinuum limit (at rn\^a = 0). We observe a good scaling be- 
havior w.rt. the lattice size, see also Fig.|5] For all couplings 
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FIG. 5: Reduced improvement term AS/N for different lattice sizes: 
9x9 (squares), 15 x 15 (triangles) and 25 x 25 (circles). Colors 
depict A = 0.8 (red), 1.0 (green), 1.2 (blue), 1.5 (magenta). 
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(red), 1.9 (blue)). 



A and m\.^a the improvement term is found to be smaller than 
14% of the total bosonic action. Depending on the coupling 
strength A, this ratio is reached sooner or later Actually, this 
represents a threshold above which the simulation fails. The 
situation is depicted in Figs. |6] and |7] At some instant, the 
improvement term blows up and settles again at a value about 
40 times the size of the bosonic action. At the same time also 
the fermion determinant grows drastically and so hinders the 
system from returning into the original (and desired) region 
of configuration space. A reason for this instability may be 
found by reconsidering the improved action 



SB = lY.hd^)^ 



W, 



(22) 
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In this form the action allows for two distinct behaviors of the 
fluctuating fields. The physically expected behavior consists 
of small fluctuations around the classical minima of the po- 
tential. Alternatively, ( |22] | allows for large fluctuations of ip 
to be compensated by large values of Wx- The latter would 
be dominated by UV contributions, and this is what we ac- 
tually observe, cf. Fig. [8] In this situation, it is definitely no 
longer possible to extract meaningful physics. Another view 
on this "broken" phase is taken in Fig. [T] While the ensem- 
ble with A = 1.4 exhibits the expected behavior at the only 
slightly larger value of A = 1.7 the simulation breaks down 
after about 5,000 configurations and for A = 1.9 the simula- 
tion is instantly found in the broken phase. 

To sum up, we have observed that the improved lattice mod- 
els may become unstable at any finite m\^a and hence any fi- 
nite lattice spacing. If and when this happens depends on sev- 
eral factors. Wilson fermions are affected in a stronger way 
while SLAC fermions remain stable for a much wider range 
of coupling constants. Apart from that, one should ensure by 
monitoring the improvement term or any other observable dis- 
cussed above explicitly that a simulation is not subject to this 
phenomenon. For the practitioner this is of course a major nui- 
sance and possibilities to avoid this matter are already under 



FIG. 7: MC history of normalized lattice mean ^i-g/m (SLAC impr., 
iV = 15 X 15, mm = 0.6, A = 1.4 (green), 1.7 (red), 1.9 (blue)). 
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FIG. 8: Mode analysis of ensembles in the physical (green, A — 1.4) 
and unphysical (red, A — 1.7) phase. Here p is the distribution func- 
tion for the modulus of the lattice momentum averaged over 25,000 
configurations (SLAC improved, A'' = 15 x 15, m\M ~ 0.6). 



investigation. Provided that one is confined to lattices smaller 
than 64 x 64 but demands the absence of finite-size effects, 
improved lattice models with Wilson fermions can be used 
for the continuum extrapolation of masses only up to A < 0.4. 
SLAC fermions can be used in the greater range of A < 1.5; 
the corresponding results will be presented further below. 



B. Setting the stage 

In Monte-Carlo simulations, importance sampling is 
only meaningful with respect to a positive measure. 
However, including dynamical fermions the measure is 
det(M) exp(— Sb). While the exponential factor is strictly 
positive (S'b is real), the positivity of the determinant cannot 
be guaranteed for an interacting theory and a possibly 
emerging sign problem has to be addressed. In order to make 
sensible comparisons with continuum calculations (which 
are most conveniently performed in an infinite spacetime) 
one furthermore must make sure that physical observables 
extracted from lattice simulations are free of finite-size 
effects. In order to check this, all simulations in this section 
are repeated in portions of fractional volume P of a fixed 
physical unit volume (with various values for / = N^a on 
a square lattice with N ~ Nt x Ns lattice points). In the 
following we consider both issues in more detail. 



1. Negative fermion determinants 

The Nicolai map in a supersymmetric theory is a change of 
bosonic variables which renders the bosonic part of the action 
Gaussian; at the same time, the Jacobian of this change of 
variables has to cancel the fermion determinant. In our model, 
this means 



(23) 



0.12 



det(A/) = det —{2{d^) + W'] 



In this light, an indefinite fermion determinant obviously cor- 
responds to a non-invertible change of variables in the contin- 
uum. 



(^ 



S,=2d(p + W'. 



(24) 



This map is only globally invertible if the superpotential is 
of degree 1 (the Nicolai map in this case has winding num- 
ber 1), i.e., for the free theory |22y. For our choice W'{ip) = 
rrup + .g(/3^ the map is not globally invertible, and there exists 
at least one point where dct(M ) vanishes iff g ^^ 0. By this 
line of argument (for the continuum formulation of the model) 
negative determinants cannot be ruled out. 

One way to cope with this in practical simulations is to use 
I det (M) I exp(— ^b) for the generation of configurations in- 
stead and to reweigh with the sign afterwards. Unfortunately, 
calculating the sign of det M is as costly as the computation 
of the whole determinant. Hence, this method becomes un- 
feasible for large lattices. A way out is to avoid reweighing 




no. 9: Probability for negative determinants (Wilson unimproved, 

iV = 14 X 14, ?nia„ = 0.43). 
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FIG. 10: Probability for negative determinants at different box sizes 
with varying lattice size (Wilson unimproved, m = 6, A = 2.0). 



within certain bounds for the parameters in which the ensu- 
ing systematic errors are negligible. Thus, we have to es- 
timate the frequency of occurrence of negative determinants 
as a function of the parameters. To obtain more reliable re- 
sults we have studied this subject with a naive inversion algo- 
rithm which computes the determinant from a LU decompo- 
sition and takes its contributions exactly into account. This is 
numerically much more involved than the standard pseudo- 
fermion algorithm, thus, this method is only applicable to 
small lattice sizes with up to 16 x 16 lattice sites. For fixed 
physical mass m it can be gleaned from Fig. |9] that configu- 
rations with a negative sign of the determinant show up only 
for A > 1.0. Furthermore, in order to understand the depen- 
dence on the lattice size and the lattice spacing we have fixed 
the coupling to A = 2.0 and run simulations on fractions P of 
a unit physical volume {I G {3/6, 4/6, 5/6, 6/6}) and differ- 
ent lattice spacings. The results displayed in Fig. [TO] clearly 
show that the problem dissolves in the continuum limit but 
becomes worse at every finite lattice spacing when the physi- 
cal volume is increased. For both figures, for each data point 



about 50,000 configurations were evaluated. Eventually, to 
estimate the impact on actual measurements we have mea- 
sured the bosonic action with m = 5, A = 2 on a 12 x 12 
lattice and obtained about 7% configurations with a nega- 
tive sign for the fermion determinant. The expectation val- 
ues considered here are (Sb) non-re„^ioh<,d = 149.94(12) and 
('^B),eweighed = 149.49(10). Hcncc even at large coupHng (far 
larger than what we target at in the next section) effects may 
be assumed to be at most of marginal relevance for actual mea- 
surements. 



2. Finite size effects 
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For these models the bare mass miatt also sets the scale for 
the overall spacetime volume. As with all lattice simulations 
we have to balance finite-size and discretization errors. If 
the lattice spacing is chosen too large, lattice artifacts may 
grow; on the other hand if, say, the Compton wavelength of 
a particle is larger than the spacetime volume the extraction 
of masses may suffer from finite-size effects. One way to 
test for the presence of such finite-size violations is to study 
the model at different spacetime volumes. Comparing the 
fermion species introduced earlier Wilson fermions may be 
expected to be most affected. Here, lattice artifacts further in- 
crease the correlation lengths so that measurements are much 
more sensitive to the finite box size. Our setup for this anal- 
ysis is as follows. At first we have simulated the improved 
lattice model using Wilson fermions at fixed coupling param- 
eters m = 15 and A = 0.3 for five different lattices with 
Nt = N, e {20,24,32,48,64} lattice points in each di- 
rection {N ~ Nt X Ns). In the following we assume that 
with this choice of coupling constants the spacetime volume 
is large enough so as to allow for a sufficiently good identifica- 
tion with the thermodynamic limit. The masses obtained from 
these simulations were extrapolated to the continuum as de- 
scribed in App.lBl This is also shown in Fig.fTTIwhere the re- 
sulting fit (and its uncertainty) is depicted with a gray shaded 
area. The next step is to decrease the volume to fractions /^ 
(with; e {9/15, 7/15, 5/15, 3/15}) of a fixed physical unit 
volume. As long as no finite-size effects are visible we expect 
to find the masses extracted at these smaller and smaller vol- 
umes to lie on top of the fit from the original lattice (of unit 
volume). Up to a volume less than half the size of the origi- 
nal one this scaling may be easily infered from Fig. [TT] which 
justifies a posteriori the correctness of our earlier assumption. 
However, since by perturbation theory the physical masses 
decrease for growing coupling (see next section), we expect 
growing Compton wavelengths and therefore stick to unit vol- 
ume (I = 1) for all further measurements so as to exclude 
finite-size effects. 



C. Weak coupling 



FIG. 11: Lattice masse.s for m = 15, A = 0.3 on box sizes I £ 
[0.2, 0.6]. We see a systematic deviation from the I = 1 result below 
1^0.5. 



the energy gap. Since unbroken supersymmetry in the con- 
tinuum predicts that bosonic and fermionic masses coincide it 
also provides a possibility to check the supersymmetric prop- 
erties of the lattice prescription. The corresponding values 
can be extrapolated from the lattice masses in the continuum 
limit. In the weak coupling it will be possible to match these 
results to predictions of perturbation theory. This provides an 
important test for the numerical results and ensures that also 
the results at intermediate coupling are reliable. 

For a description of our prescription for the boson and 
fermion mass extraction from correlators on the lattice we re- 
fer the interested reader to App. |A| With these methods we 
are aiming at a test of the lattice results against perturbation 
theory for A < 0.3. 

The reference value is given by a one-loop calculation of 
the renormalized mass 



m 1 



4A^ \ 
3737 



©(A^) 



(25) 



in the continuum valid for A ^ 1 with the bare mass m as 
used in Eq. ( fT2] i. To obtain this result one first must calcu- 
late contributions of the loop diagrams to the propagator. An 
expansion in A then yields the above result. "* 

As will be show below the fermionic masses have lower 
statistical errors than the bosonic ones. Therefore we com- 
pare only the extrapolations for fermionic masses to the per- 
turbative results. This procedure gets justified by the fact 
that bosonic and fermionic masses coincide even on a finite 
lattice for the weak coupling regime as described below in 
Sec. lIIIC3l 



An interesting observable for comparing lattice results with 
continuum physics is the mass of the lightest excited state, i.e. 



^ We will elaborate on the analytical side and the determination of the effec- 
tive potential of this theoiy in a forthcoming publication. 
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FIG. 12: The continuum extrapolation of masses for A = 0.3 for the 
improved Wilson and twisted Wilson model. Here, the SLAC result 
is given for one single lattice size. For comparison the exact results 
for the free theory are also shown. 



FIG. 13: Continuum masses for the weakly coupled regime in com- 
parison to the perturbative result. The shaded area conesponds to 
the extrapolation provided by the continuum results according to 
Eq. ^ with m = 15 and b = 1.35(13). 



1. Continuum limit 

The methods to extrapolate to the continuum given in App.lB] 
are based on the free theory with A = 0. Since we are in- 
terested in the interacting case we must first make sure that 
the continuum extrapolation of masses remains stable even for 
A = 0.3. 

To that purpose we consider the masses in the improved 
model with standard Wilson and twisted Wilson fermions at 
A = 0.3 at different lattice spacings a. In the perturbative 
coupling regime we use throughout square lattices of sizes 
Nt = N, & {20,24,32,48,64}. These correspond to lattice 
spacings of about a G [0.015625, 0.05]. A statistics of 10,000 
independent configurations puts us in a position to extrapolate 
to the continuum. 

Using these masses 7n(a) at finite lattice spacing the extrap- 
olation is shown in Fig.[T2] For comparison we also mark the 
mass for SLAC fermions at a finite lattice size N^ = N^ = 45 
(corresponding to a w 0.022). All these results indicate that 
even at A = 0.3 the continuum extrapolated masses coin- 
cide within error bounds. Even better, the masses of SLAC 
fermions at finite lattice spacing can not be distinguished from 
the continuum result. 



2. Comparison with perturbation theory 

As described above we extrapolate masses for Wilson (im- 
proved and unimproved) and twisted Wilson (improved) 
fermions for A e [0, 0.3] to the continuum values, cf. Fig. [13] 
The masses coincide within error bars although the twisted 
Wilson masses are systematically smaller. This difference has 
to be interpreted as a systematic error in the continuum extrap- 
olation for the masses but its effect is almost overshadowed by 
the statistical errors in our case. However this result indicates 
that for a reliable extrapolation at larger statistics finer lattices 



TABLE II: Continuum extrapolations of fermionic masses for Wilson 
and twisted Wilson fermions in the weak coupling regime. 



A 



Wilson unimp. 



0.02 
0.04 
0.06 
0.08 
0.10 
0.12 
0.14 
0.16 
0.18 
0.20 
0.22 
0.24 
0.26 
0.28 
0.30 



14.999(2) 
14.992(4) 
14.982(6) 
14.974(8) 
14.95(1) 
14.94(1) 
14.91(1) 
14.86(2) 
14.82(2) 
14.80(2) 
14.76(3) 
14.70(3) 
14.64(3) 
14.57(4) 
14.50(4) 



Wilson imp. 
14.997(2) 
14.993(4) 
14.999(7) 
14.963(8) 
14.96(1) 
14.91(1) 
14.87(2) 
14.87(2) 
14.85(2) 
14.79(2) 
14.72(3) 
14.73(3) 
14.60(3) 
14.60(4) 
14.45(4) 



tw. Wilson 



14.999(1) 
14.993(3) 
14.977(4) 
14.963(5) 
14.935(6) 
14.905(9) 
14.871(9) 
14.83(1) 
14.82(1) 
14.75(2) 
14.71(2) 
14.63(2) 
14.60(2) 
14.53(2) 
14.45(3) 



can be necessary to yield a better continuum limit. 

As a further test we use these results to reproduce the per- 
turbative formula 



m(A) w moWl - 



A2 



(26) 



Taken this functional form for granted, the parameters mo and 
b can be extracted from a least-square fit to the given data. 
For this fit we can use our knowledge about the free theory 
{mo ~ 15) as a fixed input or, alternatively, allow for both toq 
and b as free parameters. The corresponding results are given 
inTab.ini] 

The extrapolated results for mo confirm that the extrapola- 
tion to the free theory works reliably and that we can expect 
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TABLE III: Fit for the perturbative mass formula with 0{\^) correc- 
tions to be compared with the one-loop results. For comparison the 
one-loop result is 6 ~ 1.2990. 



derivative 


b 


mo 


Wilson improved 
Wilson unimproved 
twisted Wilson improved 


1.34(6) 
1.39(7) 
1.26(4) 


15.007(6) 
15.008(6) 
14.996(4) 


Wilson improved 
Wilson unimproved 
twisted Wilson improved 


1.37(5) 
1.42(6) 
1.25(3) 


fixed to 15 
fixed to 15 
fixed to 15 



to obtain meaningful results for b. Furthermore the results 
obtained for improved and unimproved Wilson fermions co- 
incide very well and therefore both provide the correct contin- 
uum limit. 

Additionally the results for standard Wilson and twisted 
Wilson fermions lead to compatible results when taking sys- 
tematic uncertainties of the continuum extrapolation into ac- 
count. 

As an important result of these observations, all three mod- 
els considered in the weak coupling case tend towards the 
same continuum limit for A > 0. The perturbative results can 
be recovered where the largest error bars (including possible 
systematic errors) yield b = 1.35(13) in agreement with the 
one-loop result of 6one-ioop ~ 1.2990. 



3. Signs of supersymmetry at finite lattice spacing 

Apart from all results solely based on fermions, we are pri- 
marily interested in the restoration of supersymmetry on the 
lattice. For this reason we better also check the demand from 
supersymmetry that the masses of bosonic and fermionic su- 
perpartners match. This is going to be checked by comput- 
ing bosonic and fermionic masses at couplings A = 0.2 and 
A = 0.4 with ?7i = 15 for all the models on different lattice 
sizes. 

As we have seen in the whole weak coupling regime the 
fermionic masses do not suffer from statistical noise. This 
behavior derives from the fact that the fermionic correlator 
for the free theory (A = 0) is independent of the bosonic field 
ip and is obtained by a pure matrix inversion. At small (and 
finite) A, corrections to the free propagator are of O(A^), and 
the fluctuations of if during the simulation are suppressed with 
A^; a statistics of only 10"* is needed to get reliable results. 

On the other hand the bosonic correlator even for the free 
theory is given by the correlations of the fluctuating field ip. 
Therefore a much higher statistics is necessary to sample the 
bosonic two-point function. Here, problems arise by the expo- 
nentially growing relative error of the two-point function C{t) 
with respect to t. 

Only with the use of an algorithm combining Fourier accel- 
eration with higher order integrators it was possible to simu- 
late 10^ to 10^ configurations for each parameter set [m, A) 
with an autocorrelation time of the two-point function of 



TABLE IV: For different models and lattice sizes we computed 
bosonic and fermionic masses with bare mass m — 15. 



model 



N, A 



mi 



mb,i 



mb,2 



Wilson impr. 



24 0.2 
24 0.4 



32 
32 

48 
48 



0.2 
0.4 
0.2 
0.4 



64 0.2 
64 0.4 



11.592(2) 
11.375(4) 
12.224(2) 
11.945(5) 
12.941(5) 
12.548(13) 
13.349(10) 
12.89(3) 



11.53(4 
11.39(3 
12.20(3 
11.95(3 
12.87(5 
12.47(4 
13.45(9 
12.73(9 



Wilson unimpr. 



24 0.2 
24 0.4 



32 
32 



0.2 
0.4 



48 0.2 

48 0.4 

64 0.2 

64 0.4 



11.591(2) 
11.400(4) 
12.221(2) 
11.965(5) 
12.942(5) 
12.572(14) 
13.347(7) 
12.91(2) 



11.58(2 
11.44(2 
12.20(3 
11.97(3 
12.92(6 
12.54(4 
13.45(9 
12.82(9 



11.59(4 
11.34(3; 
12.15(4 
11.88(4; 
13.02(5 
12.53(4; 
13.32(9; 
12.83(9 



11.63(3 
11.39(3; 
12.15(4 
11.87(4 
13.00(7; 
12.49(4 
13.32(9 
12.79(9 



tw. Wilson (impr.) 24 0.2 14.811(7) 14.94(11) 14.91(12) 

24 0.4 14.13(1) 14.21(9) 14.06(8 

32 0.2 14.788(6) 14.61(14) 14.94(12 

32 0.4 14.08(1) 14.39(14) 13.68(13 

48 0.2 14.789(6) 14.74(11) 14.61(11 

48 0.4 14.04(1) 14.16(16) 13.98(15 



SLAC impr. 



45 0.2 
45 0.4 



14.768(4) 
13.997(13) 



14.87(10) 
14.08(11) 



14.83(9) 
13.92(10 



SLAC unimpr. 



45 0.2 
45 0.4 



14.769(4) 
14.047(16) 



14.75(6) 
13.74(8) 



14.57(6; 
13.75(7; 



T < 2. 

The results of these numerical efforts are summarized in 
Tab. HV] They show that independently of the model even for 
A e {0.2, 0.4} bosonic and fermionic masses coiTespond to 
each other and lattice-induced supersymmetry breaking can 
not be observed. 

Finally in Figs.[T4]and[T5]the derived bosonic and fermionic 
masses are shown for the improved (and unimproved) model 
with Wilson fermions. Even these high statistics do not al- 
low for a clear cut distinction between the extrapolated con- 
tinuum masses of bosons and fermions for the improved and 
the unimproved models. This proves that even at A = 0.4 the 
improvement is not necessary even on a finite lattice. Each 
model tends towards the supersymmetric continuum limit. 



D. Intermediate coupling results 

Earlier attempts to go beyond the perturbative regime could 
not reliably determine the mass spectrum. Namely, this was 
hindered by instabilities introduced by improvement terms. 
For Wilson fermions, this renders simulations at intermedi- 
ate couplings invalid. For our analysis of coupling constants 
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FIG. 14: Bosonic and fermionic masses for the weakly coupled 
regime for the improved Wilson model. 
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FIG. 15: Bosonic and fermionic masses for the weakly coupled 
regime for the unimproved Wilson model. 



TABLE V: Fermionic masses for the intermediate coupling case. 
Twisted Wilson results are continuum extrapolations whereas the 
SLAC data is from a 45 x 45 lattice. 



A 



tw. Wilson 



SLAC unimp. 



SLAC imp. 



0.20 


14.80(2) 


14.769(4) 


14.768(4) 


0.35 


14.23(2) 






0.40 


13.99(3) 


14.05(2) 


14.00(1) 


0.45 


13.62(5) 






0.50 


13.30(6) 






0.55 


12.8(1) 






0.60 


12.2(1) 


12.81(4) 


12.44(6) 


0.65 


11.9(2) 






0.70 


10.4(5) 






0.80 




11.49(9) 


10.2(3) 


1.00 




10.2(2) 


9.4(2) 


1.20 




10.1(3) 


9.1(3) 



Apart from that, one further run was done on a 63 x 63 lattice 
with A = 0.8. Square lattices turned out to be more convenient 
with SLAC fermions and to be sufficient to clearly read off 
(within statistical errors) the masses. As for the simulations 
with twisted Wilson fermions we have determined only the 
masses from the fermionic correlators since with the statistics 
(50,000 trajectories) achieved so far the bosonic correlators 
are far too noisy to yield reliable results. 

Our results may be found in Tab. [V] and are depicted graph- 
ically in Fig. [161 From the comparison with perturbation the- 
ory first deviations are seen as soon as A > 0.4 where the 
(extrapolated) lattice results are slightly stronger curved. Also 
clear deviations between the improved and unimproved model 
using SLAC fermions become apparent for A > 0.6. It is 
worthwhile to note that the result from the improved lattice 
model is closer to the continuum limit which may be infered 



in the intermediate regime (0.3 < A < 1.2) we have there- 
fore only considered actions with twisted Wilson and SLAC 
fermions (which are anyhow expected to give better results 
at finite lattice spacing). For twisted Wilson fermions we 
have run simulations with the improved action on lattices with 
Ns e {32,40,48,56,64} lattice points in the spatial direc- 
tion. For the temporal direction we have used 1.25 • A^s lattice 
points in order to be able to assess whether contributions from 
higher excited states are really absent. At the chosen value 
of 771 = 15 in all simulations, the respective bare lattice mass 
parameter 777iatt confines the attainable coupling strengths to 
A < 0.7."* For even larger coupling strengths A only SLAC 
fermions have been found to yield sensible results. In our sim- 
ulations we used for this species both the improved and unim- 
proved lattice models on a fixed lattice size of A^ = 45 x 45. 



14.0 



12.0 




1.2 



"* For A = 0.7 we already observed that the .simulation failed on the coarsest 
lattice and had to be excluded. 



FIG. 16: Masses of the improved and unimproved model with SLAC 
fermions on a 45 x 45 lattice and continuum extrapolated results for 
twisted Wilson fermions are compared with the perturbative one-loop 
result in the continuum. 
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TABLE VI: Fermionic masses for the SLAC derivative on two dif- 
ferent lattice sizes for A = 0.8. 



N, 



improved 



unimproved 



45 
63 



10.22(26) 
10.54(15) 



11.49(9) 
10.70(19) 



from Tab. |Vll While the lattice data from the improved model 
almost coincide for both lattice spacings the data from the 
unimproved model is likely to approach the same value if in- 
creasingly finer grained lattices are used. 

Larger values of A are attainable however the numerical ef- 
fort becomes more involved and some technicalities need to 
be addressed. Once this is under control we are confident 
to investigate the strong coupling regime with the improved 
models up to A = 2.0 on the same lattice sizes. The interest- 
ing question whether the masses of superpartners still agree 
can then be satisfactorily answered. 



metric model. 

We believe that a generalization of our numerical methods 
to all supersymmetric theories without gauge fields can be ac- 
complished. In particular, the JV = 1 model in both two and 
four dimensions as well as supersymmetric non-linear sigma 
models are within reach. At least the experience gained in 
two-dimensional models suggests that SLAC and twisted Wil- 
son fermions might be good candidates for the formulation of 
four-dimensional supersymmetric lattice theories. 
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rV. CONCLUSIONS AND OUTLOOK 

In this article we have presented a detailed numerical analysis 
of the two-dimensional A/^ — (2, 2) Wess-Zumino model. Due 
to algorithmic improvements we were able to study lattice 
models at much larger lattice sizes, i.e. smaller lattice spacings 
and more importantly at stronger couplings. For a compari- 
son with analytical results from perturbation theory we have 
checked explicitly for finite-size effects and other systematic 
errors such as sign changes of the fermion determinant. Both 
were seen to be under control for the scrutinized parameter 
range. We could confirm earlier weak coupling results and for 
the first time resolve deviations from perturbation theory. All 
three kinds of fermions, Wilson, twisted Wilson, and SLAC 
fermions, approach the same continuum results. It turned out 
that lattice artifacts were largest for Wilson and smallest for 
SLAC fermions. At intermediate coupling we observed that 
the supersymmetric ally improved lattice action using Wilson 
fermions lead to unstable simulations that eventually fail to 
produce reliable results unless very large lattices are chosen. 
Simulations with SLAC fermions proved to be much more sta- 
ble; they allow for improvement terms for a wider parameter 
range. At finite lattice spacing and weak coupling no signif- 
icant differences in the measured spectrum between simula- 
tions using the improved or unimproved actions could be seen. 
It is only at larger coupling that deviations become visible, and 
the improved lattice action in fact suppresses lattice artifacts. 
It is still an open problem to go to even stronger couplings. 
Practical simulations become considerably more involved due 
to stronger fluctuations in the sign of the fermion determinant. 
Further refinements of our algorithm are already under inves- 
tigation, and we hope to report of our progress in the near 
future. Apart from that, the attainable large statistics allow 
for the determination of the (constrained) effective potential 
for this theory; this might serve as an independent check of 
the non-renormalization theorem for this particular supersym- 



APPENDIX A: DETERMINATION OF MASSES FROM 
TWO-POINT CORRELATORS 

One important observable of a quantum field theory is the en- 
ergy gap between the ground state and the first excited state. 
This energy gap corresponds to the mass of the lightest parti- 
cle in the spectrum. 

To obtain the masses in the Wess-Zumino model one has 
to consider the propagators of fermions and bosons. At van- 
ishing spatial momentum pi = 0, the free bosonic continuum 
propagator in momentum space reads 



Qboson^p^ 



1 



-Pi 



(Al) 



The real and imaginary parts ipi and ip2 of ip decouple (the 
propagatoris diagonal and even equal for (pi, ip2)- The Fourier 
transform of G^°'^°"{p) shows the well known exponential de- 
cay 



C"'°'^™(f)cxcxp(-TO|t|) 



(A2) 



where m is the above mentioned mass of the lightest parti- 
cle. (The space coordinates corresponding to pi and po are 
called X and t, respectively.) In the interacting case this quan- 
tity can be obtained on the lattice by measuring the two-point 
function. The projection onto pi = can be achieved by av- 
eraging over the spatial lattice sites, 

CrW = ^E(^"(0'0)^/3(i,x)) , (A3) 



with a, (3 labeling components of the bosonic field. 
The free fermionic continuum correlator for pi = is 



(V^a^) = G^T'°"(P0) 



*7°/3P0 



pI 



(A4) 
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Using the representation of the 7 matrices as described after 
(HJi one can read off a direct connection with the bosonic cor- 
relator using 



G' 



fermion / 



fermion / 



iPo) ■■= Gfr°"(Po) + G'™™bo) = 



2m 



-Pi 



(A5) 

As in the bosonic case on the lattice a summation over the spa- 
tial lattice sites yields the projection onto pi = 0. C''^''™'™ {t) 
defines the Fourier transform of this object. 



1. Fermion masses 

The fermionic propagator C{x) is given by 

{Mp) = (ACi(<^i: ^2)) , (A6) 

where M is the fermion matrix. The calculation of this quan- 
tity requires a high numerical effort for the inversion of large 
matrices. Fortunately in the weak-coupling limit the fermion 
matrix is approximately the same as that of the free theory and 
the statistical fluctuations are rather small. Therefore the nec- 
essary statistics to read off a reasonable fermionic correlator 
is much smaller than for bosons. 

After the fermionic correlator in position space is computed 
the masses can be determined from its long range behavior. 
Inspired by the continuum connection between fermionic and 
bosonic correlators, ( IA5b . and the behavior at large distances, 
(IA2b . one can consider 



TOeff = In 



C 



fermion 



(t) 



Cf'=™'°"(t + 1) 



(A7) 



with i in a region between zero and Nt/2. The mass can then 
be determined from the average of rricff. 

A more elaborate way is a least square fit of the fermionic 
correlator C^'=™'°"(t) with the function 



/a,m,-(i)=a-cosh(TOf(t-iV,/2)) 



(A8) 



One better not take the full range of t into account for this fit 
because it is valid only for large distances (for periodic bound- 
ary conditions, from both boundaries of the lattice). One 
should therefore constrain t to be in {1 + tskip, . . . , iVt — 1 — 
^skip}- The choice of tskip is determined by the fringe of the 
plateau in a plot of the fitting result vs. tskip- 

The differences of the different methods to determine the 
masses are illustrated in Fig. [17] One clearly observes that the 
effective masses determined according to (|A7t do not show 
a plateau from which the mass can be read off. By contrast, 
the masses obtained from a cosh fit clearly show this behavior 
at large ts^p. As mentioned above, the effective mass of the 
bosonic correlator is subject to much larger statistical errors. 
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no. 17: Bo,sonic and fermionic masse.s obtained via a cosh-fit l lASt 
and the effective mass definition l lA7t for the improved Wilson model 
with A — 0.4 on a 64 X 64 lattice. The fermionic masses with a 
statistics of about 5,000 independent configurations are much sharper 
and more reliable than the bosonic effective masses obtained from 
about lO*" independent configurations. 



the relative statistical eiTor of the correlator grows exponen- 
tially. Therefore, one must achieve a balance between this 
statistical error and the systematical errors due to the evalua- 
tion at small distances. 

We have fitted ln(C''°'™(t)) against the function A + 
\n{cosh{m\, / Ns{t — iVt/2))) to determine A and the effec- 
tive mass TTib. In order to exclude the points with the largest 
statistical and systematical errors from this fit, we have taken 
only the points in the interval ([t^kip, tst] U [Nt - tst, Nt - tskip]) 
into account, tskip is determined as in the fermionic case and 
tst such that the statistical error becomes comparably small. 

If the SLAC derivative is used an oscillatory behavior of 
nib as a function of tskip can be observed. In the bosonic 
case it is slightly smaller than the statistical error There- 
fore, it is sufficient to measure a "smeared" mass, toslac = 
0.5mb(tskip,tst) + 0.25mb(tskip + l,tst) + 0.25mb(tskip-l,fst), 
where the error of the oscillations is negligible as compared to 
the statistical one. 



APPENDIX B: CONTINUUM EXTRAPOLATION 

For the continuum extrapolation we focus on the fermionic 
masses because of their much smaller statistical error The ex- 
plicit extrapolation procedure is guided by analytic results and 
observations for the free theory. The three different discretiza- 
tions investigated in this work require different strategies for 
this procedure. 



2. Boson masses 

In order to calculate the bosonic correlators for the determina- 
tion of the masses the connected two-point function is consid- 
ered. At large distances, where the masses can be extracted. 



1. Wilson derivative 

Compared with the continuum formula, (|A5l l. the free mo- 
mentum space correlation function for the Wilson derivative 
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gets a momentum dependent mass. 



G 



fermion 



iPo) 



mua + 1 - cos(po) 



sin (po) + ("Malt + 1 - cos(po))^ 



(Bl) 



The pole of this correlator coincides with the above mentioned 
cosli-fit within the error bars. 

To extrapolate the continuum limit an expansion in powers 
of the lattice spacing is used. Exact results for the free the- 
ory were derived to check this extrapolation. In this case an 
expansion up to a linear order in a is not enough to obtain 
the known result within the high precision of the numerical 
measurements at weak coupling. Therefore we first tried to 
extended the expansion to a quadratic order which yields a 
better result; but still the error is to large for our purposes. 

The functional behavior of the masses, mf, obtained by the 
fit as a function of the lattice spacing is well approximated by 



mf{a) sa Wcont + A ■ a + B ■ a^ 



(B2) 



for all a G [0,0.05]. The deviation from this behavior is negli- 
gible with respect to the statistical errors in the weak coupling 
case. In addition the expected continuum result is achieved 
with the necessary precision. Motivated by these results this 
formula is also used in the interacting case. 



2. Twisted Wilson derivative 



A Wilson parameter of r = i/| for the twisted Wilson 
fermions in the free th eory leads to discretization errors of 
0{a^) as discussed in lUSll . For the weakly coupled regime 
(A < 0.3) we expect these errors to dominate the lattice ar- 
tifacts. Nevertheless for an intermediate coupling corrections 
of 0{a) arise. Taking this into account we extrapolate the 
masses to continuum assuming a functional behavior of 



mf{a) 



A-a + B-a^ 



(B3) 



For A > 0.3 the 0{a) terms dominate. Therefore a linear 
extrapolation is sufficient. 



3. SLAC derivative 

As we have seen in our previous investigations, lUSll . the 
SLAC-derivative shows an almost perfect behavior. That 
means the extrapolated masses coincide with their continuum 
counterparts already at finite lattice spacings. On the other 



hand we have observed an oscillatory behavior of the corre- 
lation function. This was shown to be connected with the ex- 
act reproduction of the continuum dispersion relation by the 
SLAC derivative. To handle this problem we have again stud- 
ied the free theory first. As in the bosonic case the plot of mf 
versus iskip does not show a clear plateau but rather oscillates 
around the expected continuum value, cf. Fig.fTSI 

Guided by these observation of the free theory a suitable 
averaging can lead to the extraction of the correct continuum 
results at finite lattice spacing. Starting with the ansatz 

TO(iVs, c) := comf(tsidp) + cimf(tsidp - 1) + C2mf (fskip - 2). 

(B4) 
we minimize the difference form the known continuum result 
of the free theory 



A(7V„c) = |m(iV„c)-mcont| 



(B5) 



for lattice sizes of Ns ~ Ni G {35, 37, . . . , 75} and iskip = 
[0.4A^sJ ■ A least square fit yields 

Co = 0.11791, ci = 0.47877, C2 = 0.40332 , (B6) 

leading to maxA(iVs,c) = 5.282 x 10"''. A smaller iskip 
does not change this result considerably. Using this approx- 
imation scheme the systematic error based on the oscillatory 
behavior of the SLAC derivative can be neglected compared 
to the statistical errors at least for the weak coupling case. 
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FIG. 18: Masses obtained via a cosh-iit for the free theory using the 
SLAC derivative. At larger lattices the oscillation amplitude around 
the continuum value gets smaller. 
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